Dynamical control of electron spin coherence in a quantum dot 
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We investigate the performance of dynamical decoupling methods at suppressing electron spin de- 
coherence from a low-temperature nuclear spin reservoir in a quantum dot. The controlled dynamics 
is studied through exact numerical simulation, with emphasis on realistic pulse delays and long-time 
limit. Our results show that optimal performance for this system is attained by a periodic protocol 
exploiting concatenated design, with control rates substantially slower than expected from the upper 
spectral cutoff of the bath. For a known initial electron spin state, coherence can saturate at long 
times, signaling the creation of a stable "spin-locked" decoherence-free subspace. Analytical insight 
on saturation is obtained for a simple echo protocol, in good agreement with numerical results. 
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Developing and benchmarking strategies for decoher- 
ence suppression in spin nanosystems are of vital impor- 
tance to various areas of quantum physics, from quantum 
control theory to quantum device technologies. A cen- 
tral spin 1/2 interacting with a bath of N external spins 
provides a natural testbed for detailed analysis^. This 
conceptually simple system shows a rich variety of de- 
coherence regimes, paving the way to the understanding 
of more complex scenarios, such as decoherence of many- 
spin central systems. A prominent example is an electron 
localized in a quantum dot (QD) at experimentally rele- 
vant sub-Kelvin temperatures and moderate (sub-Tcsla) 
magnetic fields, where the hyperfine coupling with a bath 
of nuclear spins is the dominant decoherence channel. Al- 
though electron spins in QDs have a broad range of po- 
tential applications in spintronics 2 - and scalable quantum 
information processing^, the coherence time T% is very 
short, T2 ^ T£, where the free induction decay (FID) 
time T 2 * ~ 10 ns in a typical GaAs QDi. While sugges- 
tive proposals exist, to increase T2 by achieving high bath 
spin polarization or bath disentanglement, or by narrow- 
ing the nuclear spin distribution 5 , methods viable in a 
wider range of physical parameters are actively sought. 

The long correlation time and distinctively non- 
Markovian behavior of the nuclear spin reservoir make 
the electron spin an ideal candidate for pulsed spin 
resonance^ and dynamical decoupling (DD) techniques 7 . 
In double-QD devices, for instance, spin singlet refocus- 
ing has been experimentally demonstrated 4 -. For a single 
QD in a large external magnetic field, Bq > 1 Tesla, 
where the nuclei simply dephase the electron spin, Hahn 
spin echoes and their Carr-Purcell-Meiboom-Gill modifi- 
cations are expected to enhance T£ by at least an order 
of magnitude in GaAs QDs&2. For weaker bias fields, 
the nuclear spin coupling induces both dephasing and 
relaxation, and the use of higher-level DD schemes has 
been invoked recentl y 10 ! 11 . However, the formal limits 
of applicability of these analyses are very restrictive and 
extremely hard to meet in practice. 

In this paper, we perform a quantitative study of the 
electron spin decoherence DD problem in regimes which 



are important for experimental DD implementations in 
QDs, yet have received little attention so far. Focusing 
on the challenging situation of zero external field, where 
dephasing and relaxation must be simultaneously elimi- 
nated, we investigate to what extent the very stringent 
formal limitations of DD methods may be relaxed. Us- 
ing exact numerical simulations, we identify promising 
DD protocols, and show that for pulse delays up to a 
factor VjV longer than naively expected from analytical 
bounds, they are still capable of extending the coherence 
time by 2-3 orders of magnitude. Special emphasis is 
devoted to the asymptotic long-time limit, where error 
accumulation is crucial and neither intuition based on 
Magnus expansion (ME) nor the quasi-static approxima- 
tion is reliable a priori. Provided that the initial electron 
spin state is known, nearly perferct coherence preserva- 
tion may be achieved for indefinitely long times. Such 
a saturation is related to the creation of a stable deco- 
herenc free subspace (DFS)i^, and may be exploited for 
stabilizing the electron spin polarization in a QD. 

Model and DD Setting. The dynamics of a single 
electron spin S coupled to a bath B of N nuclear 
spins is described by a total Hamiltonian of the form 
H = H s + H S b + H B , where H s = H S Z is the 
electron Zeeman energy in an external magnetic field 
Bo, Hsb = J2k=iAkS ■ Ifc the hyperfine contact in- 
teraction between the electron spin and the nuclei, and 
Hb = J2k>i rfci(Ifc -h — ) the intrabath dipolar cou- 
pling between nuclear spins — i 13 i 14 i 15 . S and Ifc denote 
the electron and the fc-th bath spin operators, respec- 
tively. The nuclear spin value is set to If. = l/2i£. We 
focus on the limit of zero external field Bq = and as- 
sume that the bath is initially unpolarized. The FID time 
scale is T* = (NA 2 /^)- 1 / 2 , where A = (J2 k Al/N) 1 ' 2 « 
10" 4 /ieV for typical GaAs QDs with N = 10 6 i£. Time is 
measured in units of 1/A. 

Under ideal control assumptions, DD is implemented 
by subjecting the electron spin to sequences of instan- 
taneous 7Tft rotations along appropriate control axes n, 
equally separated by the interval r. A variety of DD 
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protocols exist, based on both deterministic* 7 -*^ and 
randomize d 17 ' 18 design. In cyclic DD, the control propa- 
gator is steered through a DD group of unitary operations 
Q = {gj}, j = 0, 1, . . . , \Q\ — 1, in a predetermined order, 
as opposed to randomized DD where the future control 
path is not known in advance. Changing gi to gj requires 
the application of a DD pulse Pjj = gjg\- Thanks to the 
existence of a periodicity time scale T c = \G\t, the analy- 
sis of cyclic DD has been mostly carried out within aver- 
age Hamiltonian theory^, upper performance bounds be- 
ing determined by the dominant non-zero corrections in 
the ME for the time evolution operator. Average Hamil- 
tonian theory no longer applies to randomized DD, where 
the evolution is most directly studied in a logical frame 
that follows the applied control^ 7 -. 

Periodic DD (PDD) is the simplest non-selective cyclic 
protocol, ensuring that the unwanted evolution is re- 
moved to first order in the ME at every T n = nT c , 
n £ N, in the short T c limit. For a single spin, PDD is 
based on the irreducible Pauli group Q P = {I, X, Y, Z}1, 
which requires two-axis control sequences of the form 
C\ = CoXCqZCqXCqZ, Co denoting a free evolution 
period. Improvement over PDD may be gained by sym- 
metrized and/or concatenated design. Symmetric DD 
(SDD) guarantees that all odd terms in the ME are can- 
celled, with T c twice as long as PDD. Concatenated DDi& 
relies on a temporal recursive structure, so that at level 
£ + 1 the protocol is C e+ i = CgXCiZCgXCiZ . Here, we 
truncate the concatenation procedure at a certain level 
and repeat a periodic sequence, referred to as PCDD, 
after every A £ t (e.g., 1 — 2 leads to PCDD2). As repre- 
sentatives among stochastic protocols, we consider naive 
random DD (NRD), which corresponds to uniformly ran- 
dom pulses over Q, and symmetric random path DD 
(SRPD), where a path to traverse Q is chosen at ran- 
dom and then symmetrized as in SDD 18 . 

The use of control pulses may suit two purposes: (i) 
complete decoupling of the system from the bath, so that 
electron spin coherence is enhanced for an arbitrary ini- 
tial state; (ii) preservation of a specific initial state, in 
which case the DD sequence may be tailored accordingly. 
Two performance metrics are then appropriate. For a 
fixed initial state we use the input-output fidelity 
F(T) = Tr[p s (T)p s (0)], where p s (T) is the reduced den- 
sity operator of S at time T starting from \ifj) and trac- 
ing out the bath. For an unknown initial state, we in- 
voke minimum pure-state fidelity F m (T) — min^ F(T). 
Analytical bounds on the expected fidelity decay for var- 
ious DD protocols have been obtained for short evolu- 
tion time o 10 ' 17 ' 18 , which calls for numerical analysis in the 
long-time regime. Simulations also make it possible to ex- 
plore DD performance for values of T c beyond the strict 
convergence domain of the ME, uj c T c <C 1, where the 
highest frequency component uj c ps |^4fc | / 4 ~ NA/4. 
Let a denote the power spectrum width of the environ- 
mental coupling, 2cr re (£ fc ,4|) 1/2 = VNA 19 . We shall 
consider r — l/2cr, thus T c > 4r — y/Nu;^ 1 . To solve 
the time-dependent Schrodinger equation of the entire S 
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FIG. 1: (Color online) Minimum fidelity vs. time in the logical 
frame with r = 0.1. Hamiltonian parameters are Ho — 0, 
ro = and N = 15. For deterministic DD, data points are 
acquired at the completion of each cycle, while for NRD and 
FID this is done after every r, and for SRPD after every 8t. 
Random protocols are averaged over 10 2 control realizations. 



plus B system, we apply the Chebyshev polynomial ex- 
pansion method to the evolution operator—, and choose 
Ak > as uniformly random numbers. 

Unknown initial state. In Fig.[T]we compare F m (T) for 
the above-mentioned DD protocols. Because the charac- 
teristic time scale td for nuclear dipolar dynamics due to 
Hb is (at least) two orders of magnitude slower than the 
one due to Hsb hi typical QDs, setting Hb = is jus- 
tified for practically relevant time regimes. All schemes 
lead to substantial enhancement of the electron spin co- 
herence, PCDD2 showing the most dramatic improve- 
ment. Although, for this system both SDD and PCDD2 
remove Hsb to second order in the ME, the higher per- 
formance of PCDD2 reflects its superiority in reducing 
coherent error accumulation. The poor performance of 
NRD is expected, since its advantages over deterministic 
DD emerge only when Q is large. Contrary to the case of 
closed systems^ 8 -, SRPD does not match PCDD2 in the 
relevant parameter range, confirming the fact that irre- 
ducible DD groups and slow baths are predicted to be 
especially favorable for concatenated control 10 . 

Motivated by the above results, we proceed with a 
more in-depth analysis of the PCDD protocol. Fig. [Ha) 
compares the performance of two levels of concatenation, 
I = 2,4, for different values of r. As expected, the re- 
sults deteriorate as r increases but, interestingly, PCDD 4 
becomes worse than PCDD2. Also interesting is the ex- 
ponential fidelity decay of both protocols at long times. 
Fig. (211b) illustrates, for each value of 2ar and different N, 
the instant of time T 90 % where F m (T) for PCDD2 reaches 
90%. The results are reasonably close to each other, par- 
ticularly for larger N, supporting their applicability up 
to realistic situations with N ~ 10 6 . Lastly, we analyze 
the effect of Hb, which becomes important once the co- 
herence time is longer than 775. Let T^; be uniformly 
random numbers in [—To, To]. To avoid demanding long 
time simulations, we increase Tq manually up to values 
comparable to A^. The results are shown in Fig. [DJc), 
where a two dimensional 3x5 QD with nearest neighbor 
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FIG. 2: (Color online) (a) PCDD 2 (solid lines) and PCDD 4 
(crosses) for r = 0.3, 0.4, 0.5, top to bottom, (b) and (c): 
Tgo% vs. r for PCDD2; different bath sizes (b), intrabath 
interactions (c). Ho — in all panels, N = 15 in (a) and (c). 



intrabath coupling is considered. PCDD2 performance is 
significantly affected by a bath with fast dynamics. Al- 
though such a regime is not directly relevant to standard 
GaAs QDs, further investigation of randomized DD is 
necessary whenever Hg and H$b compete. 

Known initial state. If the electron spin is initially 
pointing along a known direction, cyclic DD protocols 
able to stabilize the input-output fidelity value for ex- 
tremely long times may be used. This is shown in Fig. [3] 
(inset), where the curves F(T) plateau after the applica- 
tion of a sufficient number n p of pulses. While asymp- 
totic saturation behavior has been reported for purely 
dephasing spin-boson models with arbitrary initial spin 
states^i, the directional dependence observed here re- 
flects the lack of a preferred direction in the error process 
generated by H$b'- & preferred direction only emerges 
through the "effective field" created by the control se- 
quence, and long-time stability depends on proper align- 
ment between such effective field and the initial state. 
In magnetic resonance language, the resulting saturation 
effect is closely related to the "pedestals" of the long- 
time magnetization signal in pulsed spin-locking exper- 
iments^. From a control standpoint, it indicates the 
dynamical generation of a stable one-dimensional DFS 
via DDi£. Consider first a selective echo protocol, say 
a single-axis PDD along the z direction, Qz = {I,Z}, 
with a corresponding (asymmetric) pulse sequence Cz = 
CqZCqZ - which we refer to as CPMG. For sufficiently 
small r, symmetrization is enforced along the z axis, as 
described by a lowest-order Hamiltonian commuting with 
Qz in the ME, and a corresponding effective field along 
zLZL; initial S'z-eigenstates are (approximate) eigenstates 
of the decoupled evolution, whereas components perpen- 
dicular to the DFS are lost in the long-time regime. 
For non-selective DD protocols based on the irreducible 
group Qp , all directions are approximately preserved for 
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FIG. 3: (Color online) Fidelity saturation for CPMG, PDD, 
SDD, and PCDD2 starting from an initial state along the half 
cycle direction. Hamiltonian parameters as in Fig. [1] The 
inset shows how the asymptotic value F sa t is reached. In the 
main panel, a number of pulses sufficient to reach saturation 
and close the cycle of each protocol is chosen, n v ~ 50. 



short times due to maximal averaging, yet long-time sta- 
bility again occurs along the direction of the dominant 
term in the ME. Notice that the latter also coincides 
with the half cycle direction of the sequence: e.g., PDD 
may be obtained from concatenation of two CPMGs, 
C\ = CxYCxY — Cy°Cx, identifying the outer y direc- 
tion as the stable one. Similarly, for SDD and PCDD2, 
the saturated components are z and y, respectively. 

Quantitative results on the dependence of the satura- 
tion value upon control parameters are given in Fig. [3] 
Interestingly, a crossing between SDD and PCDD2 occurs 
at r ~ 0.8. In particular, the CPMG sequence, which is 
not a maximal DD scheme for the Hamiltonian in ques- 
tion, leads to saturation values comparable to the other 
protocols, thus it may be useful in settings where accu- 
rate control along two axes may not be available. More- 
over, its simplicity allows for a direct analytical study 
of the saturation effect. Within the quasistatic bath ap- 
proximation (QSA)2J&, let A k = A, I = £) fc I*, and 
M = I z . After n CPMG cycles, the survival probabil- 
ity of the initial state |*(0)> = 11)® \I,M) is given by 
|(*(0)|*(2nr))| 2 = 1 - (C 2 /B 2 ) tan 2 6 cos 2 2n6, where 
C = AyJ(I - M)(I + M +1), B = A(M + 1/2), tanfl = 
d/Vl -d 2 , d = -(B/n)ahx(ftT/2), and ft 2 = B 2 + C 2 . 
Averaging over the nuclear spin bath and taking the limit 
of large n and N, F -> F sat = 1 - (1/2) / dIdMP(I, M) 

(C 2 /B 2 ) tan 2 6 with P(I,M) ~ (I/Dy/toD )e" /2 / 2D 
and D = N/4 for an unpolarized bath 19 . In the limit 
of small r, we obtain F sat = 1 - (1/16)t 2 A 2 N = 1 - 
t 2 /2T 2 * 2 . For randomly distributed A k , A 2 N i-> £ fc A 2 k . 

Fig. 2] compares the above analytical result with the 
saturation value predicted by a semiclassical approxima- 
tion, which treats the nuclear Overhauser field as an effec- 
tive random magnetic field with zero average, but finite 
variance 2 ^. The two curves superimpose, consistent with 
the fact that for large N the Overhauser field induced by 
I indeed approaches a classical field. Also shown are data 
from exact numerical simulations of a quantum spin bath 
with randomly distributed Ak- Remarkably, for short 
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FIG. 4: Fidelity saturation vs. pulse delay for CPMG from 
a known initial state. H as in Fig. [T] Circles - QSA results; 
solid line - classical random field model; dashed line with plus 
signs - exact numerical simulations. 

pulse delays the exact and QSA results are in good agree- 
ment, in spite of QSA being well known to be only valid 
for times comparable to T 2 *. Thus, (i) DD effectively ex- 
tends the region of validity of QSA; and (ii) saturation is 
entered before the QSA becomes invalid, allowing QSA 
to accurately predict F sat for short r. 

In summary, we have quantitatively characterized DD 
of an electron spin coupled to a nuclear spin bath, with 
emphasis on long-time behavior. We find that DD can 
significantly enhance the coherence time for an arbitrary 
initial state, actual performance depending on both con- 



trol and physical parameters. For a known initial state, 
the possibility of long-time saturation has been estab- 
lished numerically and analytically, which may provide a 
way for preserving the electron spin polarization without 
the need of a strong permanent magnetic field. While, 
from a practical standpoint, the estimated control time 
scales (~ 1 ns) are roughly an order of magnitude away 
from current pulsing capabilities in GaAs QDs, experi- 
mental progress is steady. In particular, single-electron 
spin rotations have been demonstrated both in gate- 
defined GaAs and self-assembled QDs 23 . In addition, 
multipulsc CPMG-DD has also been realized, not only in 
standard NMR and ESR experiments, but also in single 
solid-state centers^, which share many relevant features 
with electron spins in QDs. These advances support the 
hope that the experimental implementation of more com- 
plex protocols will be achievable in the near future. 
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